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A potential dynamics approach is developed to determine the periodic standing and traveling wave patterns 
associated with self-propelling camphor objects floating on ring-shaped water channels. Exact solutions of the 
wave patterns are derived. The bifurcation diagram describing the transition between the immobile and self- 
propelling modes of camphor objects is derived semi-analytically. The bifurcation is of a pitchfork type which is 
consistent with earlier theoretical work in which natural boundary conditions have been considered. 
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1. Introduction 


A challenge in modern-day research in biophysics and bioengineering is to create and understand bio¬ 
chemical particle systems that mimic biological cell motion. In particular, at issue is to construct particles- 
surface systems in which particles have the capability to move themselves over certain distances by con¬ 
verting chemical energy into kinetical energy |11. As pointed out in reference |11. a fundamental class 
of such man-made self-moving systems is given by objects that are floating on a medium, are driven 
by differences in the surface tension of that medium, and at the same time produce gradients of some 
substance that affects the surface tension of the medium. Irrespective of engineering applications, self- 
propelling chemical systems allow us to study principles that might be important for our understanding 
of the physics of life, in general, and the so-called active Brownian particles Ql. in particular. In fact, it 
has been shown that under certain conditions self-propelling chemical “motors” exhibit negative friction 
terms Hi that are a hallmark of active Brownian systems ana . A variety of self-propelling chem¬ 
ical motors have been studied {1,0-01 ■ In particular, the self-motion of camphor particles moving on 
water has been extensively studied, in particular, by Nakata and colleagues 115142011 . In this context, we 
would like to point out that not only single camphor particles have been considered but also the inter¬ 
action between two self-propelling camphor particles has been examined (21], the collective motion of a 
small number (about 10) of camphor particles has been investigated ]22l[2J], and the spatial and velocity 
distributions of many interacting camphor particles have been determined 

A theoretical model for the self-motion of a single, solid camphor disc floating on water has been de¬ 
veloped in terms of a Newtonian equation for the disc and a reaction-diffusion equation for the concen¬ 
tration of the surface-tension active camphor molecules on the water surface 126]. In this context, analyt¬ 
ical solutions of the concentration wave patterns under natural boundary conditions have been derived 
and it has been shown that the self-propelling mode bifurcates from an immobile mode by means of a 
pitchfork bifurcation. However, an experimentally very useful paradigm is the self-motion of a camphor 
disc on a ring channel Ilil - Eoll . 

In view of the importance of ring-shaped designs for experimental research, the present study goes 
beyond the case of natural boundary conditions. To the best of our knowledge, we will for the first time 
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present a theoretical analysis of the periodic case that can directly be applied to and compared with 
the laboratory situation. In addition, while it is plausible to assume that the chemical wave patterns of 
camphor particles on a ring exhibit a single peak, a clear explanation why this should be the case has not 
been given so far. Such an explanation will be given below. To this end, a potential dynamic perspective 
for the coupled particle-wave system will be developed. 

Explicitly, the aim of the current study is threefold. First, we will work out a potential dynamics ap¬ 
proach to wave patterns of reaction-diffusion systems in order to qualitatively discuss the possible shapes 
of the camphor concentration patterns associated with the self-motion of camphor discs floating on a ring 
channel. Second, we will derive analytical solutions for the standing and traveling concentration waves 
associated with immobile and self-propelling camphor discs, respectively, on a ring channel. Third, we 
will numerically determine the bifurcation diagram for the case of periodic boundary conditions. In do¬ 
ing so, we will show that the pitchfork bifurcation derived earlier for natural boundary conditions can 
also be found in the case of periodic boundary conditions that is frequently used in experimental re¬ 
search. 

We consider a camphor disc traveling on a ring channel filled with water. Note that the disc releases 
camphor to the water surface, which implies that in this system we need to distinguish between the cam¬ 
phor disc and the camphor concentration on the water surface. The width of the ring channel is of the disc 
diameter such that the disc can move only in one direction, along the ring. The ring has radius R and cir¬ 
cumference 2 nR. The disc position is described by the periodic variable y(f) e [0,27 tR], where t denotes 
time. Moreover, the velocity of the disc is described by v(t). The concentration u of camphor molecules 
on the water surface at time t and at a particular position x e [0, 2 tt A’| along the ring is described by the 
field variable u[x, t ) 5= 0. The dynamics of the camphor disc is given by (26] 


Tt y " 


v, 


d ay d 

— v- - ~ —— / 

df (au+ l) z ox 


ijr=y(t) 


- p v 


(1.1a) 

(1.1b) 


with a,y,p> 0. Accordingly, the disc satisfies a Newtonian equation with a force generated by the cam¬ 
phor concentration field u and a friction force proportional to the camphor disc velocity. From a mecha¬ 
nistic point of view, the camphor concentration affects the surface tension, which in turn acts as a force 
on the camphor disc (see introduction above). The effective force given as the first term on the right hand 
side of equation <l.lb) can be regarded as a gradient force of a potential, where u is the potential. That 
is, the camphor disc is driven away from regions of high camphor concentrations. The parameters y and 
a are related to certain details of the aforementioned mechanistic relationship between camphor con¬ 
centration, surface tension, and the force acting on the camphor particle, see reference [26]. Finally, in 
equation <1.11 the parameter p denotes the friction coefficient. The camphor concentration field u(x, t) 
satisfies the reaction-diffusion equation 


— u — —-u - ku + F{x- y(t)) 
dt dx z 


( 1 . 2 ) 


with k > 0. The term -ku describes the decay of the camphor concentration on the water surface due to 
dissolution and sublimation. The function F describes the increase of camphor concentration on the wa¬ 
ter surface due to the camphor disc depositing camphor molecules to the water surface. Mathematically 
speaking, F corresponds to a source term. In this context, note that we consider not too long time periods 
during which the mass loss of the camphor disc can be neglected. The source term F is defined by F - 1 
for |A| ^ r and A - x-y and F - 0 otherwise, where r > 0 is the radius of the camphor disc. Note that 
from a mathematical point of view for large r we may not think of the disc as a circular object. We may 
imagine a solid object that has the shape of a segment of the ring with segment length 2r measured along 
the curvilinear coordinate x. Then, we have r nR. In the case r — nR, the object would be a full, solid 
ring that covers the whole surface of the ring channel. 

It is important to note that equations <1.11 and <1.2) actually represent rescaled equations such that 
the variables x, y, u, and even the time variable t are given in dimensionless units. For example, the 
(real world) laboratory time h a b measured in seconds is given by the dimensionless time t occurring in 
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equations Jl.lt and J1.2) divided by a rate constant measured in 1/sec that describes the rate with which 
camphor is released from the camphor disc to the water surface |26]. 


2. Potential dynamics perspective of the chemical wave pattern associ¬ 
ated with a self-propelling camphor disc 


2.1. Potential dynamics point of view 

The objective is to study traveling wave solutions 12^ Sec. 7.3] of the form 


u(x, t ) = g(x- c t) 


( 2 . 1 ) 


of the model defined by equations Jl.lt and J1.2t . where c is the velocity of the traveling wave. Let us 
define z-x-ct as the phase coordinate of the wave. The reaction diffusion equation Jl.lt describes a 
standing or traveling wave pattern only if the source F depends on the phase coordinate z rather than 
on x and y(f). Substituting z= x-ct into A = x- y(f), we obtain A - z + ct- y{t). The requirement that 
A depends only on z implies that y(f) = c t + yo with yo e [2 nR\. This also leads to v{t) = c. The camphor 
disc moves with the same velocity as the wave pattern. Substituting equation J2.1J into equation J1.2J . we 
obtain 


d 2 d 

-r^g=-c—g-kg + F{z-y 0 ). 
OZ A oz 


Likewise, from equation (1.1) it follows that 


( 2 . 2 ) 



IdC o 

- —[flg(yo) + l] . 
ay 


(2.3) 


The position yo of the camphor disc is arbitrary. That is, if there is a wave pattern g(z, yo) with a particular 
camphor disc position yo, then yo + h has the solution g(z,yo + h) = g(z- li,yo), which shifts the wave 
pattern by h. Without loss of generality, we put yo = 0 such that 


d 2 d 

XT g^-c—g-kg + F{z) 
dz^ oz 


(2.4) 


and 


d | 

-rs 

dz lz=o 


iuc , 

— [flg(O) +1] 2 
ay 


(2.5) 


with F - 1 if z e [0, r] U [2nR-r,2nR\ and F - 0 otherwise. The wave pattern is subjected to periodic 
boundary conditions g(z) = g{z + 2nR), which implies that the relations g(0) = g{2nR] and dg(0)/dz = 
(\g[2n R) I dz hold. The traveling wave equations J2.41 and J2.5J contain the standing wave equations as 
special case. For the standing wave we put c - 0 such that 


-gX g = ~kg + F(z) 


( 2 . 6 ) 


and 



= 0. 

lz=0 


(2.7) 


These equations have a symmetric solution g(z) = g(-z). 

Equations J2.4M2.7J can be solved using the potential dynamics method for solving reaction-diffusion 
equations (see e.g. reference (3, Chap. 9]). Accordingly, g is regarded as the position of a hypothetical 
point particle that evolves in time t and moves with a particle velocity v g . In doing so, the phase coor¬ 
dinate z is replaced by the time variable t. The particle motion is subjected to a potential force with a 
potential 17(g, t ) that depends on time. Using these replacements [i.e., z — t and dg/dz — v g (t)], equa¬ 
tion <2.41 becomes 


d7 g: 


— V g = -cv g ~— V(g,t) 


d t 


dg 


( 2 . 8 ) 
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and we are looking for periodic solutions g(f) with period T — 2:tR. During the time intervals [0, r] and 
[T - r, T] we have F = 1 and the potential V gives rise to the “force” -dV/dg = kg - 1. By contrast, for 
t e (r, T - r) we have F — 0 which implies -d VIdg = kg. In total, the potential is given by 

{ -0.5k{g~ 1 lkf for t e {[0, r] u [T- r, r]} modu i us T , f? 

l -0.5 kg 2 for t otherwise. ' 


The two different potential forms described by equation (2.9) are illustrated in figure [I] and will be re¬ 
ferred to as type I and II potentials, respectively. Both potentials are inverted parabolic potentials. For 
t e {[0, r] u [T — r, T]} the potential has a peak at g = 1/A; (top panel), otherwise the potential has a peak 
at g = 0 (bottom panel). Equations (2.8) and (2.9) have to be solved under initial conditions g(f = 0) and 
v g [t = 0) that satisfy 

vM = 0) = ~ — [ag(f=0) + l] 2 , (2.10) 

aj 

see equation (2.5) . Moreover, the periodicity condition g(f + T) = g(t) implies that the boundary condi¬ 
tions g(f = 0) = g(f = T) and v g {t — 0) = v g {t- T) hold. Finally, we have the constraint g(f) 5= 0 because 
g in the original context reflects the concentration of camphor molecules. The potential dynamics sub¬ 
jected to these initial and boundary conditions is then determined by the two types of repulsive potentials 
shown in figure [I] 




Figure 1. Potentials of type I (top) and II (bottom) defined by equation (2.9) . Here k = 0.1. 

Let us discuss period solutions g(f) corresponding to standing wave patterns g(z) associated with 
an immobile camphor disc. For c — 0, equation (2.9) describes a Newtonian equation without damping 
that is solved under the initial conditions v g = 0 and g(0) = A > 0. If A > 1/k, then we have g —- oo 
for t -► oo because both potentials decay monotonously for g > 1/A;. For g(0) = 1/A;, a periodic solution 
is not possible for 0 < r < nR. However, in the limiting case r = ttR, the dynamics of the hypothetical 
particle is subjected only to the potential I force with H max = 1/A; (top panel of figure[l) such that g(0) = 
g(f) = 1/A; is the solution of equation (2.9) . That is, the wave pattern g(z) is a constant when the water 
channel is completely filled with a solid, ring-shaped camphor object. In summary, for 0 < r < nR we have 
g(0) - Ae (0,1/A;). Due to the potential I force, for t e [0,r], the particle moves “downhill” with respect 
to the type I potential V, that is, to the left towards g = 0, and reaches at time t — r a point g(r) = B. 
At t — r, the potential switches from type I to type II. Since we have v g { r) < 0, the particle continues 
to move towards g = 0. However, at this stage it moves “uphill” with respect to the type II potential 
V and, consequently, is de-accelerated. At a time point t*, the particle has zero velocity (v g {t*) = 0). 
At this instance, the particle has reached its minimal position g(F") = g m i n . Due to the impact of the 
potential II force, the particle is accelerated and starts to move “downhill” with respect to the type II 
potential V, that is, it moves to the right. Due to the symmetry of the problem at hand, the time point t* 
is half of the period: t* - 77 2- nR. When the particle moves to the right, g increases and eventually, at 
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t = 7- r, the particle reaches the position g-B again. Note that the “uphill” movement B —* g m i n and the 
“downhill” movement g m i n —>• B are described by a time-reversible Newtonian equation, which implies 
that the trajectory git) is symmetric with respect to t — Tl 2. At t — T-r we have v g > 0. In particular, 
we have u g (7- r) — -v g ir). Moreover, the potential switches from type II to type I. Due to the “initial” 
velocity v g [T— r) > 0, the particle goes “uphill” in the type I potential V, slows down, and finally reaches 
the location g = A with v g = 0 at time point t = T. In short, the periodic solution follows the sequence 

g(0) = gmax — gir) - B — g(7/2) = g min — g{T-r) = B — g(7) = g(0) 

(2.11) 

and exhibits a single maximum at g(0) and a single minimum at g(772). Moreover, the trajectory is 
symmetric with respect to t = 7/2 for f e [0, 7] (and symmetric with respect to t = 0 for t e [- 7 12, 7 12]), 
which means that the periodic standing wave pattern g(z) has a symmetry axis. 

Periodic solutions with v g (0) ± 0 related to traveling wave patterns g(z) induced by a self-propagating 
camphor disc can be discussed in a similar way. In this context, we note that a necessary condition for a 
period solution is g(0) < 1/fc again. In order to see this, let us assume that g(0) > 1 Ik holds. For c < 0, this 
implies v g (0) > 0 which implies g —>■ oo for t —► oo because both potentials I and II decay monotonously for 
g > 1/fc. For c > 0 and g(0) > 1/fc, we see that the particle has initial velocity u g (0) < 0 and moves towards 
g = 1/fc, which is the peak of the type I potential V (figure|T] top panel). If | n g (0)| is not sufficiently large, 
it will not reach the point g = 1/fc. Rather, we will have u g {t') = 0 with g(t') > 1/fc for t' e [0, 7 -] which 
implies g —► 00 to t —► 00 . That is, a necessary condition for a periodic solution starting at g(0) > 1/fc is 
that | v g (0) | is sufficiently large such that the particle passes the point g = 1 / fc during the interval [0, r]. 
However, the particle should return somehow to the initial position g(0) > 1/fc. Once the particle has 
passed the point g = 1/fc at t' e [0, r], the only way to return to the subspace g = (1/fc, 00 ) is to pass the 
point g = 1/fc at a later time point t" > t' with velocity v g (t") > 0. From git") = 1/fc and v g (t") > 0 it 
follows that g —► 00 for t —>• 00 . In summary, periodic solutions with n g (0) + 0 and g > 0 only exist for 
g(0) e (0,1/fc). By analogy to the previous discussion for the case c — f g (0) = 0, for 0 0 (camphor disc 
and wave pattern traveling to the right on the ring coordinate x ) we obtain the sequence 


g(0) = g(r) = B<A^ git*) — gmin < B — g(7- r) = C> g min 

- git**) = D = g max > A - g(7) = g(0) = >1 (2.12) 

with 0<r<t*<T-r<t**<T and v g - 0 at t* and t**. As indicated, the maximum is reached at t** e 
[7- r, 7] and does not correspond to the initial position g(0). In any case, there is a single minimum and 
a single maximum, which means that a traveling wave concentration pattern g(z) is single-peaked and 
exhibits a single minimum. The condition g{t**) > g[T) = g(0) implies for the original problem that the 
camphor disc is located to the “right” of the peak of the traveling wave concentration pattern. For c < 0, 
we conclude in a similar vein that the traveling wave pattern is single-peaked and has a single minimum. 
As anticipated above, the potential dynamic picture reveals that both the standing and traveling wave 
patterns exhibit a single peak only. However, while the standing wave patterns g(z) exhibit a symmetry 
axis, traveling wave patterns g(z) are non-symmetric. 

The next objective is to explicitly formulate the scenarios described by the sequences J2.1H and ( 12.121 
in order to obtain analytical solutions for git) and for the patterns g(z). 


2.2. Standing waves 

We will derive the next symmetric solutions of equation <2.81 with n g (0) = 0. In this case, it is suffi¬ 
cient to determine the trajectory git) in the interval [0, 7/2 = nR] because, as we will see below we can 
take advantage of the constraint v g (7/2) = 0, see equation < 12.111 . For t e [0, r], the solution git) satisfies 
equation <2.81 with c = 0, n g (0) = 0, and g(0) e (0,1/fc). The solution reads 


g(fl = 



cosh( Vlct) H —. 

fc 


(2.13) 
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In particular, we obtain 


g(r) = 


g(0)- 


coshCv^ri + i, iv( r)= Vic 

k s 


gm- 


sinh( Vkr). 


(2.14) 


For te [r,T 12] the solution git) satisfies equation 1 12.8) for c — 0 and the initial conditions g(r) and v g [r) 
given above. The solution reads 


git) 


git) + 


fgir) 


Vk 


exp 


\/k[t-r)} + - 


gir)- 


Vair) 


\/fc 


exp 


[- \/ifc(r-r)] 


= g(r) cosh f \ficit- r)| + sinh f \fk{t- r) 1 , 


which implies that 


v g {t) = \fkgir ) cosh | y/k(t-r) j + v g {r ) cosh [ Vkit- r) j 


(2.15) 


(2.16) 


The functions d2.13> — l f2TT6t involve the unknown parameter g(0). Let us use the constraint v g iT/2) = 0 
to determine g(0). Substituting t - nR with v g {t = nR) = 0 into equation <2.16) and substituting equa¬ 
tion <2.14) into equation <2.16) as well, we obtain 


Z w = — sinh f \fic{nR - r) 1 <0 
k 1 J 


(2.17) 


with £ = g(0) - Ilk and 

w = sinh( \/fcr) cosh [ \Tk{nR - r) j + cosh( \^lcr) sinh [ '/kinR - r) j > 0 . 

(2.18) 

Since the right hand side of equation <2.17) is negative and w > 0, we conclude that £ < 0. In addition, we 
see that 

1 Y f sinh [ \/1c{tiR- r)l | 

g( 0) = f+-=-<h ---(2.19) 

k /c w 

holds. From w > sinh[ \fk{7iR — r)] it follows that the factor sinh[ Vk(nR - r)] / w > 0 is smaller than 1 
such that g(0) e (0,1/fc). That is, g(0) satisfies the necessary condition for periodic solutions. 

Equations 02.13) ( f2TT5t in combination with the initial condition <2.19) describe the analytical solution 
g(f) in t e [0,TI2], As shown above, solutions of equation <2.8) with c = 0 and v g ( 0) = 0 are symmetric 
with respect to t = 0 (and t-T 12). Consequently, solutions on a full period are given by g( f) defined by 
equations <2.13) - <235t and <2.19) for t e [0, 272] and g(f) = g(-f) for t e [-772,0]. Likewise the pattern 
of the standing wave solution g(z) can be computed from equations <2.13) - <2TT5l and <2.19) for z e [0, nR] 
and g(z) = g(-z) for ze [—ttT?,0]. 

Figure [2] shows g(z) for a standing wave pattern. The solid line was obtained from equations <2.131 - 
<2.15) and <2.19) and g(z) = g(-z). The circles were obtained by solving equation <2.8) numerically in 
the interval t e [0, T12] with initial condition u g (0) = 0 by means of an Euler forward algorithm for the 
dynamics of a Newtonian particle (single time step 10 -5 ). Numerical and analytical solution methods 
showed consistent results. Although the peak of the pattern shown in figure [2] looks kinky, from equa¬ 
tions <2.13) - <2TT5l it follows that in fact the functions g(f) and g(z) are smooth functions at t = z = 0 (i.e., 
continuously differentiable). This is illustrated in the insert of figure [2] that depicts the standing wave 
pattern in a small region around z = 0. 

Let us briefly address the case r — nR. For r — nR equation <2.13) yields git) - 1 Ik. Likewise, equa¬ 
tion <2.19) reduces to g(0) - Ilk. This is consistent with the conclusion drawn above that for r - nR the 
ring channel exhibits a homogeneous concentration pattern g = 1/A;. 
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Figure 2. Standing wave solution g (solid line) computed from the analytical solution defined by equa¬ 
tions <2.131 CfTT) and (2.19) . Circles correspond to the graph g determined by a numerical solution 
method (see text). Insert shows a detail of the graph g{z) around z = 0. Parameters: k = 0.1, R = 5.0, 
r = 0.2 (as in reference 1291 ). 


2.3. Traveling waves 


Our next objective is to find solutions of equations fl2.8H2.10l for c + 0 which implies p g (0) + 0. To 
this end, we will use equations H2.8) and (2.91 to determine the trajectory g(f) and the particle velocity 
v g (t) in the full period [0, T] as function of the unknown parameters g(0), n g (0) and c. We will then use 
the three constraints given by equation A2.10) and by the two periodicity requirements g(0) = g( T) and 
v g ( 0) = i’g(T) to determine g(0), n g (0) and c. 

For t e [0, r], equation (2.8) explicitly reads d 2 g/d t 2 = -cv g + kg - 1. The solution involves the un¬ 
known parameters c, n g (0), and g(0) and reads 

g(f) = 21 i oexp{Ai t\ + A' 2,0 exp{A 2 fi + y- (2.20) 

k 


with eigenvalues 


and amplitudes 


Ai(c) — - - + 


- + k > 0, 


^2 (c) = --- 


■ + k < 0 


-Ai.o - 


g(0) 


^2,0 - 


gm 


i 

Ic 

1 

lc 


P g (0) 

A-2 

^g(Q) 

At 


1 

1 — Ai / A 2 

1 

1 — A2/A1 


( 2 . 21 ) 


( 2 . 22 ) 


In particular, we obtain 

g(r) = Ai, 0 exp{Air} +A 2i0 exp{A 2 r} + 

v g {r) = 2li, 0 Aiexp{Air} +^ 2 , 0^2 exp{A 2 r}. (2.23) 

For t e {r,T — r ), equation (2.8) explicitly reads d 2 g/df 2 = - cv g + kg. The solution involves the initial 
conditions g(r) and n g (r) listed above and reads 


g(f) = A\j exp{Ai (f - r)} + A 2 , r exp{A 2 (f- r)} 


(2.24) 


with the amplitudes 


A\ , r - 



1 

1 — Ai / A 2 



Vg[r)_ 

Ai 


1 

1 — A 2 / Ai 


(2.25) 
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In particular, we obtain 


g{T- r) = Ai, r exp{Ai(r-2r)} + A 2 , r exp{A 2 (T-2r)}, 

v g (T — r) = j4i jr Aiexp{Ai(r-2r)} + j42, r A2exp{A2(r-2r)}. (2.26) 

Finally, for te[T-r,T] equation 1 12.8) explicitly reads d 2 g/dt 2 - -cv g + kg- 1. The solution involves the 
initial conditions g(T- r ) and v g ('T - r ) (determined above) and reads 

g(f) = Ai, r _ r exp{Ai[r-(T- r)]} +A 2 ,r- r exp{A 2 [r-(r- r)]} + i 

(2.27) 


with amplitudes 


A\,T-r - 
Az.T-r = 


g [T - r) -- k 

g(r-r)-i 


VgiT-r ) 

A2 

v g iT-r) 

Ai 


1 

1 — Ai/A 2 

1 

1 — A 2 / Ai 


In particular, we obtain 

g(T) = A\j~ r expfAi r} + A 2 ,r-r exp{A 2 r} + , 

fi 

v g (T) = ^i,r-rAiexp{Air} +A 2 ,r-rA 2 exp{A 2 r}. 

/2 

As mentioned above, we put g(0) = g(T) and v g (0) = v g (T) such that 

g(0) = /i(g(0), v g (0), c) , n g (0) = / 2 (g(0), v g (0),c) . 


(2.28) 


(2.29) 


(2.30) 


The functions f\ and f 2 occurring in equation ( 12.30) are defined by the expressions on the right hand 
sides of equation <2.29) and by the amplitudes and initial conditions listed in table [l] These amplitudes 
and initial conditions are explicitly defined by equations <2.22) . <2.23) . <2.25) . <2.26) and <2.28) with the 
eigenvalues Ai and A 2 given by equation <2.21) . 

The two relations in equation <2.30) together with the constraint <2.10) provide three equations to 
determine the parameters g(0), i/ g (0), and c. Once g(0), u g (0), and c have been determined, the trajectory 
g(t) can be computed from equations <2.20) . <2.24) and <2.27) . In doing so, the analytical solution for the 
traveling wave pattern g(z) can be found. It can be useful to eliminate c with the help of equation <2.10) 
like c — -ay n g (0)/{/x[ag(0) + l] 2 }. Equation <2.30) then becomes 

g(0) = /n(g(0), n g (0)) , v g (0) = h 2 (g(0), i/ g (0)) (2.31) 


Table 1. Amplitude and “initial conditions” involved in the definition of the functions f\ and / 2 occurring 
in equations <2.291 and <2.301 . 


Amplitudes 


Initial conditions 

A\,T- r [g{t-r), v g {T-r),c] 
A 2 ,T-rlg(t- r), Vg{T- r),c\ 


g{T -r)[A u ,A 2 , r ,c] 
Vg(T- r)[Ai : r,A 2j ,c\ 

Ai, r [g(r), v g [r),c] 
A 2 ,rlg{r),v g {r),c] 


g{r)[A lfit A 2 fi,c] 

Vg[r)[Ai' 0 ,A 2t o,c] 

Ai l0 [g(0),i/ g (0),c] 

^ 2 , 0 [g(0), f g (0), c] 
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Figure 3. Traveling wave solution g (solid line) computed from the analytical solution defined by equa¬ 
tions <2.201 . <2.241 and <2.271 . Circles were obtained by the same numerical solution method as used in 
figure [2] Parameters: a = 0.5, y = 2.0, k = 0.1, R = 5.0, r = 0.2, p = 0.1 (as in reference ||2^|). Estimated 
parameters obtained both from MATLAB fsolve and the u>i(t) and dynamical system method: 
g(0) = 0.3151, v g {0) = 0.1516, c = 1.1315. 


with 


hi = fi[-,c = -ay Vg(0)I {p[ag(0) + l] 2 }), (2.32) 

hi = -ay v g (0)/{p[ag(0) + l] 2 }). (2.33) 

We numerically solved equation 1 12. 311 for a given set of parameters a, y, k, R, r, p in order to ob¬ 
tain g(0) and fg(0). Subsequently, we calculated c from c — -ay n g (0)/{p[ag(0) + l] 2 }. We solved equa¬ 
tion <2.311 in two different ways. First, we used the MATLAB solver (f solve) for coupled nonlinear func¬ 
tions. Second, we numerically solved a dynamical system that exhibits a fixed point consistent with equa¬ 
tion <2.311 . More precisely, the two-variable dynamical system defined by dw\/t)/dt = w\ - hiiwi, wi) 
and dm 2 (r)/df = hiiwi, w 2 ) - u >2 was numerically solved (Euler forward with single time step 10 -5 ) for 
initial conditions wi(0) = 0.5/k and wii 0) = -0.1 (to obtain solutions with c > 0). For both methods, the 
same numerical values for g(0), v g { 0), and c were obtained up to 4 digits after the decimal point. Subse¬ 
quently, the traveling wave pattern g(z) was computed from equations <2.201 . <2.241 and <2.271 . The result 
is shown in figure[3] A detail of the pattern for values of z close to zero is shown in the insert. We see that 
the camphor disc at z = 0 is located to the “right” of the peak of the concentration wave (i.e., the peak 
position is to the “left” of z = 0). 

Finally, we compute the bifurcation diagram for the same set of parameters a, y, k, R, r as used 
to calculate the wave pattern g(z) in figure [3] However, p was not fixed. Rather, the parameter p was 
considered as control parameter and was varied in the range p = [0,0.2]. The traveling wave velocity c as 
function of p was determined using the two aforementioned methods (MATLAB f solve and dynamical 
system uq(f), W 2 {t)). The propagation velocities c thus obtained are shown in figureO The bifurcation 
diagram indicates that for periodic solutions g(z) there is a pitchfork bifurcation at a critical value p c 
such that at p = p c , the standing wave solution with c - 0 becomes unstable and two stable traveling 
wave solutions emerge with either c > 0 or c < 0. 


3. Discussion 

A potential dynamics approach was used to determine the shape of chemical concentration patterns 
on a ring-shaped water channel associated with a camphor disc floating on the water surface of the 
channel either in an immobile or self-propelling fashion. The potential dynamics approach allowed us 
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Figure 4. Pitchfork bifurcation describing the transition of a camphor disc floating on a water ring from 
an immobile mode to a self-motion mode. Two methods were used: MATLAB fsolve (solid line) and 
dynamical system (f) (circles). See text for details. 


to qualitatively determine the shape of the wave patterns and to show that the concentration patterns 
are bounded from above. In rescaled dimensionless units, the maximal camphor concentration of stand¬ 
ing or traveling wave patterns cannot exceed the inverse of the decay constant k. With the help of the 
potential dynamics approach, analytical expressions for standing and traveling wave patterns were de¬ 
rived. Finally, by means of the analytical expression for traveling wave patterns, the bifurcation diagram 
describing the transition from standing to traveling waves (immobility to self-motion) when the friction 
parameter /i is decreased was determined numerically. 

The stable branches of the bifurcation diagram for values of /i smaller than the critical value /i c but 
close to /i c (i.e., /i = /i c - e with e > 0 small) may be described like /i c - /i = Lc 2 + 0(c 4 ), where L > 0 
holds and c corresponds both to the camphor disc velocity and the traveling wave velocity. Taking the 
immobile-disc and standing-wave solution c - 0 into account, this relationship between /i and c becomes 
(/i c - jJ) c = Lc 3 + 0(c 5 ) and describes a supercritical pitchfork bifurcation. This relationship was analyt¬ 
ically derived in a previous study assuming natural boundary conditions [2jJ. Consequently, our consid¬ 
erations indicate that a supercritical pitchfork bifurcation can also be observed in the case of periodic 
boundary conditions. 

In the aforementioned study it was also shown that under natural boundary conditions and when the 
radius of the camphor disc becomes large, the velocity c and the control parameter fj. can be related like 
(ju c - fi) c = L\c 3 + L 2 C 5 + 0(c 7 ) with L\ < 0 and L 2 > 0. When the solution c = 0 is ignored, the simplified, 
truncated relation /i c - /i = L\ c 2 + L 2 C 4 describes the branches of self-propelling camphor discs associated 
with traveling wave concentration patterns. It has the shape of a W when considering e — fi c - /i as a 
function of c. Consequently, this second relationship describes a subcritical pitchfork bifurcation when 
considering c as a function of e — - /i. The subcritical bifurcation predicts that for appropriately chosen 

values of p, there are two different propagation velocities 0 0 (and likewise two velocities c < 0). We 
were not able to find a counterpart of this observation for the case of periodic boundary conditions. 
That is, the two numerical methods described in section [2~3l to determine the bifurcation diagram did not 
indicate the existence of a W shaped relationship between ^t c - /i and c. In this context, it is important to 
note that in reference [26j] it was shown that the range of values fi for which two different positive (or 
negative) velocities c exist is relatively small. Therefore, the issue of a subcritical pitchfork bifurcation 
for solutions of the reaction-diffusion equation dl.2l subjected to periodic boundary conditions may be 
addressed by more sophisticated analytical or numerical methods that go beyond the scope of the present 
study. 
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TOMHi p03B'fl3KM XBI/IJlb XilVliHHOI KOHLjeHTpaLJM 
caMoypyxoM/itoBa/ibHMX i<aM<()opHMX Hacn/iHOK, U40 pyxatoTbcn 
no Ki/ibL^io: HOBa nepcneKTi/iBa noTeHi^ia/ibHoV Ai/maMiKi/i 

T./],. OpaHK 

Biflfli/i ncuxonorii, YmBepcuTeT m. KoHHeKimyr, CT 06269, CLLIA 

P 03 Bl/ 1 HyT 0 MeTOfl nOTeHpia/lbHOI flMHaMiKW TOTO, 11406 BMBHHTH KapTHHH nepiOflWMHUX CTOflHHX i 6i>KyHWX 
XBn/ib, noB'a 3 aHwx i 3 caMoypyxoM^K)Ba/ibHHMM KaMcfiopm/iMn o6'eKTaMH, axi pyxaioTbCfl Ha xi/ibi4eno,qi6HMX 
BOfl^Hwx xaHa/iax. OTpwMaHO TOHHi po 3 B'fl 3 KH XBW/ibBOBHX xapTWH. OTpMMaHO HaniBaHa/iiTWHHO fliarpaMy 6i- 
cJjypKapiY, axa onwcye nepexifl Mix HepyxoMora i caMoypyxoM/uoBa/ibHoio MOflaMM xaMcfiopHi/ix o6'exTiB. Bitfiyp- 
xaL(m e BW/ionofli 6 Ha, mo y 3 roflxyeTbC 5 i 3 nonepeflHboto TeopeTHHHOto po 6 oTOto, b nxiPi po 3 r/iflHyro npi/ipoflHi 
rpaHWMHi yMOBH. 

K/iiOHOBi cnoBa: XBiwi xiMiHHo'i KOHU,eHTpau,ii, caMopyx, Bi^jionofli6Ha 6icpypi<apifi 
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